Exact Solutions of Domain Wall and Spiral Ground States in Hubbard Models 



Makoto Homma and Chigak Itoi 
Department of Physics, Nihon University, Kanda, Surugadai, Chiyoda, Tokyo, Japan 

We construct a set of exact ground states with a localized ferromagnetic domain wall and an 
extended spiral structure in a deformed flat-band Hubbard model. In the case of quarter filling, we 
show the uniqueness of the ground state with a fixed magnetization. We discuss a more realistic 
situation given by a band-bending perturbation, which can stabilize these curious structures. We 
study a conduction electron scattered by the domain wall and the spiral spins. 
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Domain wall and spiral structures in ferromagnetic sys- 
tems are interesting structures. Generally a domain wall 
is localized stably between two ferromagnetic domains. 
In this case, the ferromagnetic order is preserved within 
one domain, and the translational symmetry is sponta- 
neously broken. Domain structures are believed to ap- 
pear universally in ferromagnetic systems with an energy 
gap and with a finite correlation length. On the other 
hand, spiral structures appear in special situations. Such 
structures are extended over the entire space, and some- 
times destroy ferromagnetic order. In this letter, we con- 
struct a set of exact ground states in a class of Hubbard- 
like models with coexisting domain wall and spiral struc- 
tures, which have never been found in other models. 
Some remarkable results for ferromagnetic ground states 
have been obtained in a class of Hubbard models recently. 
Mielke and Tasaki showed independently that the ground 
state gives saturated ferromagnetism in a many-electron 
model on a lattice with special properties, which is called 
the flat-band Hubbard model [jj, |2j . Tasaki also proved 
the stability of the ferromagnetism against a perturba- 
tion which bends the electron band . Tanaka and Ueda 
have shown this stability for a two-dimensional model in 
Mielke's class |^|. If domain wall structures are univer- 
sal, then a domain wall solution within this framework 
should give us some important physical insight into the- 
ories of many-electron systems. Here, we construct such 
a solution. We deform a flat-band Hubbard model by in- 
troducing a complex anisotropy parameter q. The SU(2) 
spin rotation symmetry in the original flat-band model is 
reduced to U(l) in our deformed model. This anisotropy 
| q | 7^ 1 leads to a localized domain wall with a finite 
width and a complex q leads to an extended spiral state. 
We study the stability of the domain wall ground states 
against a band-bending perturbation using a variational 
argument. We prove that for \q\ ^ 1, the energy ex- 
pectation value of a state with a domain wall centered 
near the origin becomes lower than the eigenvalue of the 
saturated ferromagnetic eigenstate, unlike the SU(2) in- 
variant model. We discuss similarities of the wall solution 
and differences between the domain wall solution of our 
model and domain walls in pure quantum spin systems. 



Alcaraz, Salinas and Wreszinski constructed a set of ex- 
act ground states with two domains in the XXZ model 
with a critical boundary field in arbitrary dimensions for 
an arbitrary spin [f|. They showed that the degener- 
acy of the ground states corresponding to the location of 
a domain wall center is identical to that of the ground 
states in the SU(2) invariant model. In their solution, the 
domain wall is localized at an arbitrary surface with a fi- 
nite width depending on the Ising anisotropy parameter 
(q + q~ 1 )/2 > 1. We will see that the domain wall ground 
state in our electron model has the same degeneracy as 
that in the XXZ model and the same localization prop- 
erty in a certain parameter regime. On the other hand, 
for complex q, our model differs from the XXZ quantum 
spin model which has no spiral ground state. Finally, we 
study from a microscopic viewpoint the scattering of a 
conduction electron by the domain wall and spiral spins. 

Here, we consider a one-dimensional lattice with a site 
index x = —L — 1, • • • , L + 1, where x is an integer and 
L is an odd integer. Electron operators on the lattice 
satisfy the anticommutation relation 

{Cxa, 4t} = SarSxy, <7, T =f, 4 • 

We define the following electron operators for even site x 
with a complex number q and a real number A > 

dx-it = -(<7 1/4 )*Cx-2T + Ac x _i T - (q~ 1/4 )*c xh 
a x -n = -{q~ 1/4 )*c x -2i + Ac^-ij. - {q 1/4 )*c xl , 
dx] = q~ 1/4 c x -i^ + Ac^t + q 1/4 c x+lh 
dxi = q 1/4 c x -n + Ac^x + q~ 1/4 c x+ n, (1) 
n x1 = c 1t C:e T' = c tl Cx i' Ux = nx 1 + nx -i- ( 2 ) 

The operators a xa and <F yT are always anticommuting 

{a xa ,dl T } = 0, x = odd, y = even, a, r =|, J. . (3) 
The Hamiltonian is written in terms of the above opera- 



2 



tors as 



H — H nop + Hi n t (4) 
H hop = t (4f^T + 4:i d xl) = E ^cLcV ( 5 ) 



x,y,a 



(6) 



a zero-energy eigenvalue is a ground state. An arbitrary 
state created only by is a zero-energy states because 
of its anticommutativity @ with d x i and no double oc- 
cupancy, and thus this polarized state 



|allup)= JJ 



vac), 



(8) 



with a repulsive coupling constant U > and a hopping 
parameter t > 0. The anticommutativity and the absence 
of a x<T in H nop give a band flatness of the electrons created 
by a X(T . Here we use an open-boundary condition under 
which there are no degrees of freedom at the edge sites 
x = —L — 2 and x = —L+2, namely c_l-2<t = = cl + 2<t- 
The hopping amplitudes on a unit cell defined by H nop 
are depicted in Fig. ^ 
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FIG. 1: Arrows and circles depict hopping amplitudes and 
on-site potentials, respectively, for electrons with up(a) and 
down(b) spins on a unit cell. For complex q = | g| e 10 , the hop- 
ping amplitude opposite to the arrow is the complex conjugate 
of that along the arrow. 

Let us discuss the symmetry of this Hamiltonian. 
The first important symmetry is U(l) symmetry. The 
Hamiltonian commutes with the total magnetization 
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Pauli matrices 
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S&\ where we define = J2x S ^ and S ^ = 

a}cxp, where crW (j = 1,2,3) are the 
In the limit q — > 1, this symmetry is en- 
hanced to the SU(2) symmetry. In this limit, this model 
becomes the original flat-band Hubbard model given by 
Tasaki 0, ■ The second important symmetry is defined 
by a product of a parity P and a spin rotation 



n = n^ 1 =Pcxp (mS^ 
nc^fjll C—X —(j ^Y\.c xo .TL - - c_ x 



(7) 



Note the following transformation of the total magneti- 
zation nS^^n = — S( 3 \ An energy eigenstate with the 
total magnetization M is transformed by n into another 
eigenstate with the total magnetization —M, which be- 
longs to the same energy eigenvalue. 

Now we construct ground states in the quarter-filled 
case, where the uniqueness of the ground state can be 
shown. Since each term in the Hamiltonian is positive 
semi-definite H nop > and -ff; n t > 0, an eigenstate with 



is a ground state. To create ground states with other 
magnetization, first we should take into account a condi- 
tion of no double occupancy at odd sites. A state created 
by a xa . under this condition is written in the following 
summation over all spin configurations with an arbitrary 
fixed magnetization ^ x a x = M 



\M)= E 



1p(<T-L, ■ 



n 



vac). 



-M 



odd 



The condition of no double occupancy on the state \M) 
at even sites yields 

^(••- ,T,l,---) = ^(--- ,!,?,•••)• 

This relation implies the uniqueness of the ground state 
with a fixed total magnetization, since two arbitrary spin 
configurations can be related by the successive exchanges 
of two nearest neighbor spins. Therefore the degeneracy 
of those ground states is exactly the same as that in the 
SU(2) symmetric model, as in the domain wall solution 
in quantum spin systems |5|. 

To explore the nature of the ground state, we write it 
in a more explicit way. The following superposition over 
the states with different magnetizations can be a ground 
state 



E^ 

M 



Mi 



M >= 



11(4 



x1 +zqza xl 



f )|vac) 



(9) 



a;— odd 



where z is an arbitrary complex number. This state can 
be regarded as a generating function for a ground state 
with an arbitrary magnetization. The expectation values 
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FIG. 2: Spin one point function on odd sites x evaluated from 
the numerical solution of the recursion relation for L = 1000, 
q = 1.05, A = 3 and z = 1.05" 
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of the spin operators at site x in this ground state are where r = (e + Ve 2 — 4)/2. The convergence of the quo- 



9(1) \ = zg 2 + (z<p)* 
1 / 2 l + |zV| 



(2)\ _ ( n x) m* 2 - N 2 )* 

2i l + \z 2 q x \ 
{n x )l-\z 2 q x \ 



a(3)\ _ \"x/ \* H I 

where the expectation value of an operator O is defined 
by (<PL= (^l k)/( z l z >- As discussed in the XXZ mod- 
els 0, S @ j the two domains are distinguished by the 
sign of the local order parameter (S^). The domain wall 
center is defined by the zeros of (S^) which is located 
at x = — 21og| g | \z\. The function \{n x ) — |(>si 3 ^)| decays 
exponentially from the center, if (n x ) ^ 0. The domain 
wall width 1/ log |g| is defined by its decay length. A pro- 
file of the domain wall in a typical example is depicted 
in Fig [3 For the complex q = \q\e , one can see the 
spiral structure with a pitch angle 9. The vector (S x ) is 
rotated with angle 6x around the third spin axis depend- 
ing on the site x. Note that this spiral structure of the 
ground state does not exist in the XXZ model. Though 
the complex anisotropy parameter q = e lS is possible in 
the XXZ Hamiltonian, the spiral state is no longer the 
ground state and the corresponding model is described 
in the Tomonaga-Luttinger liquid without ferromagnetic 
order. The translational symmetry in the infinite volume 
limit is broken by the domain wall or the spiral structure 
for finite log \z\. Both symmetries generated by and 
IT are broken spontaneously as well. 

Now, we evaluate the function (n x ) to see the practical 
form of (S x ) given by eq. (|10|l . First we consider a simple 
case of the large A limit. This limit implies that there is 
no overlapping of the electrons at each even site, thus 

(n x ) = 1 - 0(A~ 2 ), x = odd, (n x ) = 0(A~ 2 ), x = even. 

^ (11) 

The localization property of the ground state in this case 
is the same as in the XXZ model. To treat a general case, 
we introduce the following normalization function 

X 

A{x 0l x,0 = \\ nKr + ^rM^lvac)!! 2 , (12) 

y=xa 

where we employ z as a real number and use a parameter- 
ization z 2 = \q\- c for \q\ > 1. Note A(-L,L,() = (z\z). 
This normalization function obeys the recursion relation 

A(x , x, =e (1 + \q x - c \)A(x ,x - 2, C) 

-(l + l^-^D^o^-^c), (13) 

where e = A 2 + \q\ 1 ^ 2 + |g| -1 ^ 2 - We can extract the main 
x dependent part out of A(xq, x, () as follows 

X 

A(xo,x,()=B(x 07 x,Or*-^ ±1 (1 + 1^1), (14) 



tient £>(— L, a;, £) in the infinite volume limit is proved 
rigorously on the basis of the recursion relation for 
B(—L, x, C) derived from the recursion relation (|13|l 10]. 
The one point function is represented in terms of the 
normalization function 

_ A 2 B(-L,x-2,QB(x + 2,L,Q 
[Ux) ~ r B(-L,L,Q 

for odd x. The expression for even x is also easily ob- 
tained as a closed but complicated form. We have bounds 
I < (n>x) < u with some positive constants I and u less 
than 1. Therefore, these bounds and eq. jjHfll guarantee 
the unique domain wall center at x = £. The numerical 
solution of the recursion relation for B(—L,x,Cj is useful 
for seeing the profile of the domain wall practically. The 
spin one point function for {S x ^) is depicted in FigEl 

Now, we study the stability of the domain wall and the 
spiral structure of the ground state against a band bend- 
ing perturbation in a variational argument. To bend the 
lower flat band, we consider the following perturbation 



(15) 



with a positive small constant s. As has been rigorously 
proved for q = 1 by Tasaki |3|, the all-spin- up state © 
is preserved as a ground state for sufficiently small s/t 
and s/U. Here, we consider the problem also in this 
limit. The all-spin- up state (JBJl is an eigenstate of the 
Hamiltonian H = Hhop + -Hint + -f^hop still 

ff|all up) = £J(all up) | all up), 

where £7 (all up) = —{L + l)se. If \q\ ^ 1, other do- 
main wall states © are not eigenstates, and the energy 
expectation values are lower than this energy eigenvalue 
E(al\ up). First, we consider the case of \q\ ^ 1. In the 
large t and large U limit, the energy expectation value is 
written in terms of the normalization functions 



E(0 



{Am 

(z\z) 



E(allup)-s ^ g(x,C), (16) 



— L 



where x is summed over odd integers and 
(\q\h-\q\-h) 2 B(-L,x-i,0 



r 2 (\q^\ + \q-^\)B(~L,LX) 
2eB(x + 2, L, () (\q\ \ + \q\-^) 2 B(x + 4, L, Q 



(17) 



rQq 1 



|g 2 I) 



Let us compare the energy expectation values of two so- 
lutions with different domain wall centers, namely,^ and 
C+2. The analysis of the recursion relation in ^(| enables 
us to estimate B(—L, x, £), then the energy difference of 
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two nearest neighbor domain walls with a finite distance 
£ = O(l) from the origin is 

\E(( + 2)-E(0\=O(L\q\- L S ). 

The energy cost by the finite shift of the domain wall 
center might vanish in the infinite volume limit. If this is 
true, the domain wall solution with an arbitrary center 
is stable against the perturbation . The degeneracy 
corresponding to the location of the domain wall center 
may be preserved in the infinite volume limit. On the 
other hand the energy difference between the all-spin-up 
state and the state with a domain wall at C = — L + 1 is 

£(all up) - E{-L + 1) = O(s). 

The energy eigenvalue of the all-spin-up state is larger 
than the domain wall state with order 1. Since the en- 
ergy of the true ground state is lower than the expecta- 
tion value of the trial state, the all-spin-up state is no 
longer a ground state unlike the SU(2) invariant model. 
Therefore, in this variational argument, we conjecture 
that the localized domain wall structure is stable against 
this band-bending perturbation. On the other hand, for 
\q\ = 1, all the states \z) are still exact eigcnstates and 
are still degenerate to the eigenvalue .©(all up). There- 
fore, the stability of the ground states can be proved in 
the same way as that given by Tasaki 3] . 

Next, we consider the behavior of one electron added 
to the quarter filled ground state \z) with a domain wall 
or a spiral structure, 

According to a variational principle, we find the trial 
wave function ip xa to optimize the following variational 
functional 

"rTp^^ T ~TTra^^J (na;> - (19) 

In our model, the scattering of conduction electrons by 
the domain wall or spiral spins of the lower band has 
been given already as a well-defined problem. By opti- 
mizing this effective energy function under a normaliza- 
tion constraint (ip\ip) = 1, we obtain a two-component 



Schrodinger equation. At large A or for \q\ = 1, the (n x ) 
becomes periodic for the lattice unit cells. In particular, 
at large A, eq. Ijlljl implies that a localized spin sits at an 
odd site and the conduction electron hops mainly on even 
sites at large U. The normalization constraint becomes 
local in this limit and the hopping term gives a conven- 
tional kinetic term in continuum approximation, and the 
equation for real q > 1 becomes that for a conduction 
electron on the Kondo lattice coupled to the XXZ model 
obtained by Yamanaka and Koma The model on 

the Kondo lattice is realized under special conditions in 
our model. 



In this letter, we construct a set of exact ground states 
with a ferromagnetic domain wall and a spiral structure 
in a deformed flat-band Hubbard model. Some results 
obtained here can be extended to higher dimensional 
models, as discussed in the original flat-band Hubbard 
model Q. After the construction of the domain wall 
ground state in the XXZ model JyL interesting properties 
of excitations were discovered US IE E3] I n particular, 
singular low-lying excitations localized near the domain 
wall were discussed in arbitrary dimensions 0,0. We can 
expect similar properties of excitations above the domain 
wall in our electron model as well. 
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